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Abstract. We use WR124 (WN8h) and its associated 
nebula Ml-67, to test theoretical non-LTE models for 
Wolf-Rayet (WR) stars. Lyman continuum ionizing flux 
distributions derived from a stellar analysis of WR124, 
are compared with nebular properties via photo-ionization 
modelling. Our study demonstrates the significant role 
that line blanketing plays in affecting the Lyman ionizing 
energy distribution of WR stars, of particular relevance to 
the study of H n regions containing young stellar popula- 
tions. 

We confirm previous results that non-line blanketed 
WR energy distributions fail to explain the observed neb- 
ular properties of Ml-67, such that the predicted ionizing 
spectrum is too hard. A line blanketed analysis of WR124 
is carried out using the method of Hillier & Miller (1998), 
with stellar properties in accord with previous results, ex- 
cept that the inclusion of clumping in the stellar wind 
reduces its wind performance factor to only ~2. The ion- 
izing spectrum of the line blanketed model is much softer 
than for a comparable temperature unblanketed case, such 
that negligible flux is emitted with energy above the He I 
A504 edge. Photo-ionization modelling, incorporating the 
observed radial density distribution for Ml-67 reveals ex- 
cellent agreement with the observed nebular electron tem- 
perature, ionization balance and line strengths. An al- 
ternative stellar model of WR124 is calculated, follow- 
ing the technique of de Koter et al. (1997), augmented 
to include line blanketing following Schmutz et al. (1991). 
Good consistency is reached regarding the stellar proper- 
ties of WR124, but agreement with the nebular properties 
of Ml-67 is somewhat poorer than for the Hillier & Miller 
code. 

Key words: stars: Wolf-Rayet - starsdundamental pa- 
rameters - stars:individualWR124 - ISM: HII regions - 
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1. Introduction 

The Lyman ionizing energy distributions of hot, massive 
stars are important in the study of young starburst regions 
and galaxies, via population synthesis codes (e.g. Leitherer 
et al. 1999). However, the interstellar medium (ISM) con- 
spires to prevent this energy from reaching the Earth's 
atmosphere. Interstellar dust severely degrades transmit- 
ted fluxes at ultraviolet and optical wavelengths, while the 
Lyman continuum of hot stars is completely absorbed by 
intervening hydrogen atoms for all but a few nearby B 
stars (e.g. Cassinelli et al. 1995). Therefore, indirect mea- 
surements of the extreme UV flux are necessary 

Until recently, the standard method of obtaining the 
ionizing distributions of hot stars was to rely on line blan- 
keted, plane parallel, LTE model predictions from Ku- 
rucz (1991). However, the need to consider non-LTE ef- 
fects, spherical geometry and line blanketing effects has 
recently led to the development of complex model atmo- 
spheres (e.g. Hubeny & Lanz 1995; de Koter et al. 1997; 
Pauldrach et al. 1998), with Lyman continuum ionizing 
flux distributions obtained via synthesis of the accessible 
UV and optical stellar spectrum. However, these different 
codes predict quite different ionizing spectra for early-type 
stars of identical temperatures. While attempts at testing- 
theoretical stellar O-type models via observations of stel- 
lar clusters have been made (Stasihska & Schaerer 1998; 
Oey & Kennicutt 1998), definitive results require studies 
of individual stars and their associated H II regions. 

In the case of Wolf-Rayet (WR) stars, evolutionary 
synthesis models for young starbursts known to contain 
these stars, known as 'WR galaxies', still rely on un- 
blanketed, non-LTE model atmospheres of Schmutz et al. 
(1992). What effect does line blanketing have on the pre- 
dicted Lyman continua of WR stars? Recent theoretical 
calculations by Crowther (1999) suggest that the effect can 
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Table 1. Summary of basic stellar and nebular properties 
of WR124/M1-67. 



Quantity Measurement 



Reference 



Massey (1984) 
Massey (1984) 
van der Hucht et al. (1985) 
Crowther et al. (1995b) 
Crawford & Barlow (1991) 
Smith et al. (1996) 

Grosdidier et al. (1998) 
N e , T e 1050 cm" 3 , 6200K Esteban et al. (1991) 

Fn a 2.84X10" 11 ergcm" 2 s _1 Grosdidier et al. (1998) 



m v 11.58 mag 

(b - v) 0.81 mag 
mK 7.71 mag 
Eb~v 1.18 mag 
Distance 4-5 kpc 
Sp. TypeWN8h 

R 55-60 arcsec 



be substantial, such that unblanketed WR models overes- 
timate the hardness of their ionizing spectra, which is sup- 
ported by observations of extra-galactic giant H n regions 
by Bresolin et al. (1999). 

In principle, Hn regions associated with their central 
star are ideal tracers of the Lyman ionizing flux distri- 
butions. However, suitable nebulae are rare, since their 
properties are often poorly known. The principal stud- 
ies for WR stars are those of Rosa & Mathis (1990) and 
Esteban et al. (1993) who combined WR model fluxes 
from Schmutz et al. (1989, 1992) with observed proper- 
ties of WR ring nebulae, to investigate the properties of 
the central stars. Esteban et al. varied the temperatures 
of unblanketed WR models until agreement was reached 
between the observed nebular properties and those pre- 
dicted by photo-ionization modelling. In general, compar- 
isons with (independent) stellar analyses of the central 
stars was found to be reasonable, except that lower tem- 
peratures were required from the photo-ionization models 
for late-type WN (WNL) stars, especially the Galactic 
WN8 star WR124 and its associated Hn region, Ml-67, 
for which no agreement was achieved. 

In this work, we shall depart from the technique of Es- 
teban et al. in that stellar properties of WR124 are fixed 
from a spectroscopic study, which are then tested against 
the nebular properties of Ml-67 via photo-ionization mod- 
elling using CLOUDY (Ferland 1996). Sirianni et al. (1998) 
showed Ml-67 to be very clumpy. Hubble Space Tele- 
scope (HST) imaging of Ml-67 by Grosdidier et al. (1998) 
revealed an astonishingly complex nebula, allowing mea- 
surements of the nebular radial density distribution and 
integrated Ha flux. Our work therefore represents the first 
attempt to test the reliability of predicted WR Lyman 
continuum distributions using both unblanketed (Hillier 
1987) and blanketed models (following Hillier & Miller 
1998, 1999) plus robust nebular properties. In our study 
we compare these predictions with the code of de Koter et 
al. (1997), augmented to allow for line blanketing follow- 
ing Schmutz (1994, 1997). In this way, we will be able to 
quantitatively assess the effects of different codes and line 



blanketing techniques on the ionizing energy produced by 
WR stars. 

In Sect. 2 we will discuss observations of WR124 and 
-Ml-67. In Sect. 3 the model atmospheres codes are intro- 
-duced and discussed, with stellar results compared in Sect. 
4. The photo-ionization modelling technique is discussed 
and results are presented in Sect. 5. Finally, conclusions 
are drawn in Sect. 6. 

2. Observations of WR124 and its associated 
nebula, Ml 67 

We discuss the observed properties of WR124 and Ml- 
67 in this section, and provide a summary in Table [l]. 
Unlike most massive WR nebulae, which are difficult to 
identify, Ml-67 was so striking that it was included in the 
Bertola (1964) Planetary Nebula catalogue (PK50+3 1) 
with WR124 as its central star. Cohen & Barlow (1975) in- 
stead proposed a massive WN origin, which was supported 
by Solf & Carsenty (1982) and Esteban et al. (1991), who 
found nebular abundances to be consistent with processed 
stellar ejecta, indicative of a massive central star. Con- 
firmation of a non-PN origin was made by Crawford & 
Barlow (1991) who estimated a distance of 4-5 kpc from 
interstellar Nai D2 observations. (In order for Ml-67 to 
be a PN, its distance should be — 460pc according to van 
der Hucht et al. 1985). 

2.1. WR124 

WR124 (alias 209 BAC, He 2-427) is a Galactic, heavily 
reddened WN8h star that is well known as the ionizing 
star of Ml-67, with a high heliocentric recession velocity 
of -200 kms- 1 (Merrill 1936). Crowther et al. (1995b) 
presented optical spectroscopy of WR124, including com- 
parisons with other late-type WN (WNL) stars. 

For the current analysis, spectrophotometry obtained 
in 1991 September at the 2.5m Isaac Newton Telescope 
(INT), covering AA4400-7300 (resolution 2-3A) are used, 
taken from Crowther et al. (1995b). Observations ob- 
tained at the Deutsch-Spanisches Astronomisches Zen- 
trum (DSAZ) 2.2m observatory in June 1991 provide an 
additional blue dataset covering AA3400-4400 (~2A reso- 
lution) , details about which can be obtained from Hamann 
et al. (1995). 

New near-IR spectroscopy of WR124 was obtained in 
1998 July at the 3.8m U.K. Infrared Telescope (UKIRT) 
with the cooled grating spectrograph CGS4, the long 
(300mm) camera, a 256x256 InSb array and a 401/mm 
grating. First order observations provided data covering 
1.66-2. 30^m (R~800) and 2.98-3.62/xm (R~1300) with 
a slit width of 0.6" (1 pixel). This data set was bias- 
corrected, flat-fielded, extracted and sky-subtracted us- 
ing CGS4DR (Daly & Beard 1992). Subsequent reductions 
and analysis were carried out using FIGARO (Meyerdierks 
1993) and DIPSO (Howarth et al. 1995). In order to re- 
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Fig. 1. Comparison between the spectrophotometry of 
WR124 (open circles: IUE; solid line: INT) and theoreti- 
cal CMFGEN models, reddened by £?b-v=1-0, 1.15, 1.3 and 
1.45 mag (dotted lines) according to Cardelli et al. (1989). 
For clarity, comparisons are offset vertically by 0, +2, +4 
and +6 units. 



move atmospheric features, the observations were divided 
by a standard star (whose spectral features were artifi- 
cially removed) observed at around the same time and 
similar air mass. In regions of low atmospheric transmis- 
sion at UKIRT the reliability of line shape and strength 
must be treated with caution (e.g. Paschen a). 

The reddening to Ml-67 is extremely high, with pre- 
vious estimates of E-Q-y in the range 0.90 (Esteban et 
al. 1991) to 1.50 (Solf & Carsenty 1982). Analysis of 
the stellar spectrum by Crowther et al. (1995b) indicated 
i/B-v = l-18 mag, which implied a distance of ^6kpc, as- 
suming a typical M v =— 6.0 mag for WN8 stars. Here we 
attempt an improved stellar reddening determination by 
including International Ultraviolet Explorer (IUE) newly 
extracted spectra (INES, Rodriguez-Pascual et al. 1998) of 
low resolution (LORES) datasets. The quality of LORES 
INES datasets are preferred to NEWSIPS (Schartel & 



Skillen 1998). Four large aperture IUE datasets of WR124 
were obtained between 1978-1986, including 3 long wave- 
length (LWR) datasets, totalling 270 minutes exposure 
time, plus 1 short wavelength (SWP) dataset of 212 min- 
utes. 

Pre-empting the results from the next section, we com- 
pare spectrophotometry with reddened, line-blanketed 
model atmosphere distributions obtained with the CMF- 
GEN code (see Sect. 3T) in Fig. |l[ indicating 2?b-v~1-3 
mag (Ay=4.0), 0.1-0.2 mag higher than previous results 
based solely on optical stellar datasets. A comparison with 
alternative ISA-wind (see Sect. |3~2| ) model fluxes indicates 
an equivalent reddening. 

Our measurement is in reasonable agreement with the 
nebular study of Ml-67 by Chu & Treffers (1981) who 
obtained i?B-v=l-23 from a comparison of its 1.4 GHz 
and Ha flux. In contrast, Esteban et al. (1991) obtained an 
significantly lower value of £"b-v— 0-90 mag from nebular 
Balmer line strengths. 

Since the reddening towards WR124 is so high, we pre- 
fer to select the distance from the K-band absolute magni- 
tude, Mk, where the interstellar reddening is much lower. 
Using WN8-9 stars at known distance in our Galaxy or 
the Large Magellanic Cloud (LMC) as calibrators (see 
Morris et al. 1999 and references therein) we obtain 
Mk=— 6.2±0.3 mag. Using the observed K-band magni- 
tude of WR124 from van der Hucht et al. (1985) plus A K = 
0.12 xA v = 0.49 (Cardelli et al. 1989) we find a distance 
modulus of 13.44 (i.e. 4.9 kpc), in accord with Crawford 
& Barlow (1991). Consequently, we adopt a distance of 

5 kpc to WR124, i.e. Af K =~6.25 mag. Note that the re- 
sulting V-band absolute magnitude of My— —6.1 mag is 
in excellent agreement with LMC WN8-9 stars (Crowther 

6 Smith 1997). 

2.2. Ml-67 

Ml-67, discovered by Minkowski (1946), has a clumpy, 
ejecta-type, morphology (Chu 1981), with an angular di- 
ameter of 90-120 arcsec (Solf & Carsenty 1982; Grosdidier 
et al. 1998). High spatial resolution spectroscopy of Ml-67 
was obtained by Esteban et al. (1991), revealing nebular 
conditions of A^^IO 3 and T e ~6,200K. Relative to normal 
Hn regions, oxygen is deficient, with nitrogen enhanced, 
such that N/O («3). Their results suggest that O has been 
processed into N mainly through the ON cycle, and con- 
firm that the Ml-67 nebula is composed of stellar ejectra. 
Esteban et al. estimated an upper limit of 0.2 for the ratio 
between the ISM mass swept-up by Ml-67 and the nebu- 
lar mass, suggesting that relatively little mixing with the 
local interstellar gas has occurred. 

Fabry-Perot interferometry (Chu & Treffers 1981) and 
long-slit spectroscopy (Sirianni et al. 1998) indicate a com- 
plicated velocity field, with two different motions: a spher- 
ical, hollow shell, expanding at 46 kms -1 , plus a 'bipolar' 
outflow with a larger velocity of 88 kms -1 . Sirianni et 
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al. (1998) interpreted this structure as the result of two 
subsequent outbursts. Grosdidier et al. (1998) used the 
very high spatial quality of the Wide Field and Planetary 
Camera 2 (WFPC2) aboard HST to reveal filamentary, 
chaotic substructures, with a number of 'bullets' in the 
nebula, which seem to concentrate within two conical re- 
gions in agreement with the previous results of Sirianni et 
al. (1998). Globally, Grosdidier et al. (1998) found that the 
de-projected radial density distribution of Ml-67 is well 
described by a power law, N e (r) ~ r -°- 8 j with a cutoff at 
r >50". 



3. Stellar atmosphere codes 

In this section we introduce and utilise the codes that 
are used to carry out spectral synthesis of WR124. CMF- 
GEN (Hillier 1987, 1990; Hillier & Miller 1998, 1999) solves 
the transfer equation in the co-moving frame, subject to 
statistical and radiative equilibrium, assuming an expand- 
ing, spherically-symmetric, homogeneous or clumped, at- 
mosphere. Populations and ionization structure are con- 
sistent with the radiation field. 

To take into account the (microscopic and macro- 
scopic) velocity structure expected in a real WR star, a 
uniform Doppler line width of Vdop=50 kms -1 is assumed. 
Since this is substantially greater than the (species depen- 
dent) thermal Doppler width, we refer to it as a 'broad- 
line' approximation. All lines within some multiple of Voop 
can interact within the Sobolev resonance zone. Because 
of the large line width used, the interaction between two 
lines may be overestimated in the broad-line approxima- 
tion. 

A simplifying 'super level' approach is used for indi- 
vidual levels (Anderson 1989), particularly for iron-group 
elements. In this approach, several levels of similar ener- 
gies and properties are treated as a single 'super level', 
with only the populations of the super level included in 
the solution of the rate equations. The population of an in- 
dividual atomic level in the full model atom is determined 
by assuming that it has the same departure coefficient as 
the corresponding super level to which it belongs. 

Unfortunately, individual line blanketed co-moving 
frame calculations are generally computationally demand- 
ing, despite the use of super levels, so we first need to 
confirm that individual case studies, such as this, are in 
accord with observation. A computationally quick method 
is to (i) solve the transfer problem in the Sobolev approx- 
imation rather than the co-moving frame, which the code 
ISA-wind does (de Koter et al. 1993, 1997), and (ii) con- 
sider line blanketing via Monte Carlo sampling following 
Schmutz (1994, 1997) allowing the opacity of a huge num- 
ber of lines to be considered. How consistent are results 
obtained with these different methods? 



3.1. CMFGEN 

The model calculations are based on the iterative tech- 
nique of Hillier (1987, 1990) named CMFGEN. Allowance 
is made for line blanketing and clumping following the for- 
mulation of Hillier & Miller (1998) to which the reader is 
referred for specific details. The model atom contains hy- 
drogen, helium, carbon, nitrogen, silicon and iron as shown 
in Table 0. Details of each ion are included, such that for 
hydrogen, n=l to 30 full levels are considered, which are 
grouped into 10 super levels. Test calculations indicated 
that ~5000 iron transitions with gf < 10~ 4 have neg- 
ligible influence on the emergent spectrum, but required 
considerable computational effort, so these have been ex- 
cluded. Within a super-level, representing a group of full 
levels, populations are assumed to be in LTE. Elemental 
abundances other than hydrogen and helium are fixed at 
cosmic (Si=0.1% by mass, Fe=0.2% by mass) or approxi- 
mately CNO-equilibrium (C=0.03%, N=l to 2% by mass) 
values for a solar metallicity environment. 

The stellar radius ) is defined as the inner boundary 
of the model atmosphere and is located at Rosseland opti- 
cal depth of ~20 with the stellar temperature (T») defined 
by the usual Stefan-Boltzmann relation. Similarly, the ef- 
fective temperature (T e g) relates to the radius (-R2/3) a t 
which the Rosseland optical depth equals 2/3. 

3.2. ISA-wind 

The improved Sobolev approximation code (iSA-wind) is 
described in detail by de Koter et al. (1993, 1997). The 
principal differences with CMFGEN relate to: (a) the treat- 
ment of the line radiation transfer; (b) the wind elec- 
tron temperature, which assumes a grey, LTE atmosphere 
rather than radiative equilibrium, where the temperature 
is fixed at a pre-assigned value in the outer regions, such 
that lines originating in the outer wind (e.g. He 1 A10830) 
may be incorrectly predicted; (c) the velocity structure 
deep in the atmosphere; (d) the specific atomic model 
treated, as listed in Table ||; (e) the neglect of clumping. 

ISA-wind solves the line radiation transfer using an 
improved version of the Sobolev approximation. In co- 
moving frame codes, line photons with a finite (given) 
width interact with the continuum over a range of depth 
and frequency points, while in the Sobolev approxima- 
tion the continuum opacity and source function are con- 
sidered constant within the line resonance volume. In the 
improved version of the Sobolev code, ISA-wind, absorp- 
tion of line radiation by the continuum is taken into ac- 
count within the resonance volume. The Sobolev approx- 
imation should be ideal for the large velocity gradients of 
WR outflows. It introduces great simplification in the rate 
equations, and makes the overall iteration process about 
10 times faster than a co-moving frame code, for the same 
number or ionic levels and grid points. 
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Table 2. Summary of model atom used in CMFGEN and ISA-wind radiative transfer calculations, including full levels 
(Np), super levels (TVs) and total number of transitions (Arr ans ), following Hillicr & Miller (1998) and dc Koter et al. 
(1997). Note that transitions considered in CMFGEN include those between individual terms in multiplcts. The Monte 
Carlo calculation used to derive line blanketing accounts for transitions of most metals, including Al, Si and Fe. 
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Turning to line blanketing, an iterative technique in- 
cluding the Monte Carlo method of Schmutz et al. (1991) 
and Schmutz (1994, 1997) is used. The method allows the 
computation of intensity-weighed effective opacity factors, 
which account for the presence of tens of thousands of 
spectral lines, dominated by Fe and Ni. Based on ISA-wind 
atmosphere calculations, the Monte Carlo (MC) code de- 
termines the line blanketing factors. An iterative proce- 
dure is used, such that blanketing factors are used by the 
non-LTE code to calculate a new atmosphere, which in 
turn is used to calculate new blanketing factors. A few 
iterations are generally sufficient. This is due to the fact 
that the scattering and absorption factors are not very 
sensitive to the specific model parameters. 

The MC method deals with the radiative transfer in 
the correct manner, except that the ionization and exci- 
tation equilibrium of metal species is approximate. This 
dictates which lines are efficient at capturing photons for 
each point in the atmosphere. The ionization structure 
of H and He in the MC code is derived from ISA-wind. 
For other species, an representative ionization tempera- 
ture is used, obtained from the principal ionization stages 
for each metal that is considered in ISA-wind. The MC 



code calculated ionization stratification can therefore be 
compared with that of ISA-wind for consistency. 

Since line-line interactions are not accounted for in the 
MC calculation, explicit correction factors, named 'photon 
loss' factors, need to be considered to allow for interaction 
between important resonance lines (Hen A303, Hci A584) 
and nearby metal lines (Schmutz 1997). We have calcu- 
lated photon loss factors from the opacities of He n A303 
and lines at nearby wavelengths, and find that line-line 
interactions are negligible for the stellar temperatures of 
WNL stars. (Pasquali et al. (1997) previously adopted a 
factor of 1(T 4 at Hen A303 for WNL stars.) Additionally, 
we have determined that line-line interaction between He I 
A584 and adjacent metal lines would not affect the ioniza- 
tion balance or derived stellar parameters for WR124. 

4. Stellar analysis 

Unblanketed non-LTE spectroscopic analyses of WR124 
by Hamann et al. (1993) and Crowther et al. 
(1995a) revealed T* «30-34kK, log(L/L )~5.3-5.4, 
^f~10- 4 - 2 M Q yr- 1 , and ^-700 kms" 1 . The latter study 
included hydrogen, helium, carbon and nitrogen, deriving 
H/He^0.6 by number. Hamann & Koesterke (1998a) as- 
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Fig. 2. Comparison between the rectified optical (INT) and near-IR (UKIRT) spectrum of WR124 (solid lines) with 
the synthetic spectrum produced by CMFGEN, accounting for line blanketing and clumping. 



sumed this hydrogen content, and derived equivalent stel- 
lar parameters based on an analysis of the nitrogen spec- 
trum. In contrast, Esteban et al. (1993) used properties of 
Ml-67 to conclude that T*<30kK, although a suitable un- 
blanketed model was not found from the grid of Schmutz 
et al. (1992). The best agreement was obtained with the 
only line blanketed models available at that time, namely 
those for the WN9 star R84=Brcy 18 (T*=28.5kK) from 
Schmutz et al. (1991). A spectral comparison of WN8 and 
WN9 stars is presented by Crowther & Smith (1997). 



4-1. Analysis technique 

Our approach follows that of Crowther et al. (1995a), such 
that diagnostic optical lines of Hci (A5876), Hen (A4686) 
and H I (H/3+He n A4859) are chosen to derive the stellar 
temperature, mass-loss rate, luminosity and hydrogen con- 
tent. In the absence of high resolution UV observations, a 
wind velocity of 710 kins -1 is obtained from optical Hci 
P Cygni line profiles. A (3=\ velocity law is adopted, since 
this reproduces the optical, near-IR and mid-IR spectra of 
a similar WN8 star, WR147 (Morris et al. 1999). 
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Fig. 3. Comparison between the rectified optical (INT) spectrum of WR124 (solid lines) with the synthetic spectrum 
produced by ISA-wind, accounting for line blanketing. Spectral comparisons in the red and near-IR are comparable to 
CMFGEN, except for those lines discusssed in the text. 



4-2. Results 

In Fig. ||, we compare spectroscopic observations of 
WR124 with the synthetic spectrum obtained from CMF- 
GEN, allowing for both line blanketing and clumping. The 
level of agreement between optical observations and the- 
ory in Fig. H is excellent, except that the Hei absorption 
components are often predicted to be too strong. 

We also include comparisons with H and K band 
UKIRT observations of WR124 in Fig. § A gain, excel- 
lent agreement is reached. Note that UKIRT observations 
of Hen 3.09^im are also well reproduced. The case of 
He I 2.058/xm deserves special comment. Crowther et al. 
(1995a) discussed the difficulties in predicting emission at 
He I 2.058/xm for the WN8 star HD 86161, in that the pre- 
dicted emission was far too weak. The strength of 2.058/im 
(and He I A5016) is very sensitive to the optical depth 
of the He I A584 transition. Crowther et al. (1995a) at- 
tributed this failure to the neglect of line blanketing in 
the extreme-UV. Our results support this claim for the 
case of WR124, the predicted 2.058/xm emission now is in 
excellent agreement with observation. 

Crowther et al. (1995a) obtained a similar level of con- 
sistency in their study, based on an identical optical data 
set, except that electron scattering wings were strongly 
overestimated because of their assumption of homogeneity 
(their Fig. 10). Allowing for clumping takes away this dis- 
crepancy, such that electron scattering wings are well re- 
produced with an assumed volume filling factor of /=10%, 
with a corresponding decrease in the required mass-loss 
rate by a factor of ^3. As discussed by Hillier & Miller 
(1999) and Hamann & Koesterke (1998b), a unique de- 
termination of the filling factor is not possible, such that 
an uncertainty of ±50% remains on the precise mass-loss 
rate. 

The stellar parameters for this model are compared 
with those from the unblanketed, homogeneous analysis 
of Crowther et al. (1995a) in Table |. We find that the 




T*/kK 

Fig. 4. Contour lines of equal rectified line peak intensity 
for He i A5876 (4.3, solid) and Hen A4686 (3.55, dashed), 
in the stellar T»-log M plane. The thick contours repre- 
sent the observed values, while the thin contours show the 
±10% error margin. The rhomboids mark the parameters 
for which models were calculated. 



inclusion of line blanketing has only a minor effect on the 
derived parameters for WNL stars. Nevertheless, also con- 
sidering the effects of clumping on the mass-loss rate, the 
wind performance factor, Mw 00 /(L H «/c), is reduced from 
10 to -2! 

Turning to ISA-wind, we present a contour plot for the 
dependence of our diagnostic He i-ii lines on temperature 
and mass-loss for WR124. We prefer to use a compari- 
son in terms of peak intensity rather than line equivalent 
width since the former approach was followed for CMF- 
GEN. Our figure allows us to quantify uncertainties in pa- 
rameters, namely ±1000K in temperature and ±0.1 dex 
in mass-loss rate. Since these represent formal uncertain- 
ties with ISA-wind, how do the derived stellar parameters 
compare with results from CMFGEN? Table |3| reveals that 
good agreement is obtained between the two codes in de- 
rived stellar parameters, including the hydrogen content. 
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Table 3. A comparison of stellar parameters obtained for 
WR124 using CMFGEN (unblanketed and blanketed) and 
ISA- wind (blanketed). The unblanketed results are taken 
from Crowther et al. (1995b), based on identical obser- 
vations. In all standard velocity law with (3=1 is 
assumed. 



Parameters 


CMFGEN 


ISA-wind 


Filling factor, / 


100% 


10% 


100% 


Blanketed? 


no 


yes 


yes 


T, (kK) 


33.5 


32.7 


32.8 


R* (Re) 


17.1 


18.0 


16.2 


log(L./L Q ) 


5.52 


5.53 


5.44 


log(M/v7) Moyr" 1 


-4.1 


-4.2 


-4.4 


v-o (kms 1 ) 


710 


710 


710 


H/He 


0.6 


0.7 


0.8 


N/He 


0.004 


0.006 


0.003 


log Q (s" 1 ) 


49.05 


49.05 


49.00 


log Qi (s" 1 ) 


48.19 


43.88 


47.91 



This comparison demonstrates the validity of the Sobolev 
approximation for WR analyses, and represents an impor- 
tant result from the present study. 

In Fig. |^, we compare identical optical observations of 
WR124 with the synthetic spectrum obtained from the 
ISA- wind contour plot, allowing for line blanketing. The 
quality of line fits is comparable to CMFGEN, except that 
(i) the region around A4100 is more severely underesti- 
mated, due to the neglect of Siiv A4088-4116, (ii) He I 
A5015 is too weak; (iii) Hen A4200, A4542 and A5412 are 
in better agreement than CMFGEN. In the near-IR, Hei 
2.058/im emission is predicted to be somewhat too strong, 
as is He I 2.112/im emission, with Br7 too weak. 

4-3. Comparison of ionizing flux distributions 

Our three solutions indicate very similar stellar properties 
for WR124 and are equally successful at reproducing its 
optical spectrum. How do the Lyman ionizing flux distri- 
butions of these models compare? 

The unblanketed CMFGEN ionizing energy distribution 
is shown in Fig. ^. It is relatively flat, with a small jump 
beyond the He I edge at 1.81 Ryd (24.6 eV or 504A), and 
a non-ncglig ible flux up to 3.2 Ryd (286A). Fig. | reveals 
that the blanketed CMFGEN distribution is much softer 
than the unblanketed solution, showing a large deficit for 
energies >1.3Ryd (i.e. A <700A), such that a very low flux 
is emitted at energies beyond the Hei edge. In addition, 
this distribution shows a moderate excess between the Ly- 
man edge and 1.25 Ryd relative to the unblanketed case, 
which explains why its bolometric luminosity is essentially 
identical. 

ISA-wind shows an intermediate energy distribution, 
with a depression around 1.7 Ryd (540A), but a relatively 
hard flux beyond the Hei edge, up to energies of 2.9 Ryd 




.5 1.0 1.5 2.0 2.5 3.0 3.5 




-2 - 



.5 1.0 1.5 2.0 2.5 3.0 3.5 

Energy (Ryd) 

Fig. 5. (upper panel:) Comparison between the Lyman 
ionizing distribution of unblanketed CMFGEN (dotted 
lines) and blanketed CMFGEN (solid lines) for WR124, re- 
vealing large differences; (Lower panel) as above, except 
for unblanketed CMFGEN (dotted lines) and blanketed ISA- 
wind (solid lines). Note that these diagrams are designed 
to emphasise the extreme UV spectral region, important 
for the nebular ionization, and not the overall energy dis- 
tribution. 

(310A). Possible explanations for the differences in ioniz- 
ing fluxes between the two blanketed cases are discussed 
in Sect. |t[ 

5. Photo-ionization modelling 

From the usual diagnostic diagram relating Ha/[Sli] to 
[S n] A6717/A6731 (see Sabaddin et al. 1977) we find that 
Ml-67 falls in the photo-ionization dominated region. 
Therefore we use the general purpose photo-ionization 
code CLOUDY (Ferland 1996) for our analysis. The ion- 
ising radiation fields presented in the previous section are 
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expected to yield different predicted nebular properties, 
when input into CLOUDY, which can be quantitatively 
compared with spectroscopic nebular observations taken 
from Esteban et al. (1991). 

5.1. Description of the calculations 

We constructed photo-ionization models using CLOUDY 
(v90.04) as described in Ferland (1996) and Ferland et 
al. (1998). Comparisons with other photo-ionization codes 
are provided by Ferland et al. (1995). The nebula are rep- 
resented by a sphere of variable gas density, n, and filling 
factor, e, with a central cavity that is ionized and heated 
solely by the UV radiation of a single central star. Nebular 
fluxes are predicted, given input abundances, flux distri- 
butions and physical parameters. We use a volume filling 
factor of e=0.05 as estimated by Solf & Carsenty (1982). 
The de-projected radial density profile obtained by Gros- 
didier et al. (1998) was used, together with an outer (in- 
ner) nebula boundary at 50" (8"), which corresponds to 
1.25 pc (0.2 pc) using a distance of 5 kpc. 

Grosdidier et al. (1998) provided a measurement of 
the integrated nebular Ha flux, after correction for [Nil] 
contamination. We obtain an integrated de-reddened Ha 
flux using the reddening of £'b-v=1-3 mag, determined 
in Sect. |2.1| Tfm™\— k /iiwin-io 



namely /(Ha)=5.45x 10 _iu ergcm^s -1 . In 
addition, long slit nebular line intensities of Esteban et al. 
(1991) (their position A) are de-reddened according to this 
value of Eb-v- Although the majority of previous stel- 
lar and nebular reddening determinations are in reason- 
able agreement with our value, Esteban et al. (1991) ob- 
tained a substantially lower value of £b-v=0.90 mag for 
Ml-67 using Ha/H/3. Consequently, the Ha/H/3 ratio, de- 
reddened according to our stellar reddening of £b-v = 1-3 
is far from the usual Case B value. 

In order to assess the impact of alternative reddening 
determinations on the nebular properties, we have also 
calculated the integrated Ha flux, following the reddening 
determination of Esteban et al. (1991), i.e. £b-v=0.9, so 
that /(Ha)=2.35x 10 -10 ergcm _2 s _1 . For this case, long 521 CMFGEN 



Abundances are assumed to be typical of Galactic H II 
regions, except for N, O and S, which are adapted from 
Esteban et al. (1991). For each reddening, electron densi- 
ties and temperatures were obtained from the de-reddened 
[S II] A6716/A6731 and [N 11] A5755/A6584 ratios using the 
ratio program (Adams & Howarth, priv. comm.). For 
E B -v=Q-9 we obtained iV e =10 29 cm~ 3 and T e =6,000K, 
while N e =10 30 cm^ 3 and T e =6,200K for £ B -v=l-3. Sub- 
sequently, ionic abundances were obtained by solving the 
equations of statistical equilibrium using the equib pro- 
gram (Adams & Howarth, priv. comm.), and are given in 
Table |j for each reddening. The usual elemental diagnos- 
tics were chosen, [Nil] A6584, [On] A3727, [Sn] A6731 
and [S ill] A9069, with no ionization correction factors 
for nitrogen (since 2+ <C + ). Also, we have adopted 
He/H~0.2, since helium is expected to be moderately en- 
riched. Solar abundance ISM dust grains were assumed, 
although the precise choice is not critical. 

5.2. Results from photo-ionization modelling 

Esteban et al. (1993) identified a large discrepancy be- 
tween the nebular fluxes of Ml-67 as determined by un- 
blanketed model atmospheres (Schmutz et al. 1992) and 
that observed. In particular, they were unable to repro- 
duce the observed T e and S + /S 2+ ratio. Closest agree- 
ment was achieved for the line blanketed model of R84 
(WN9) from Schmutz et al. (1991). We now investigate 
whether this discrepancy remains, based on improved stel- 
lar models, including those for which line blanketing is in- 
corporated. Comparisons between predicted and observed 
nebular properties are presented in Table |J, including re- 
sults from both the reddening of Esteban et al. (1991), or 
the present value, to illustrate the possible range of line 
strengths and plasma conditions. Lamers et al. (1999) dis- 
cuss various aspects of the reliability of generally derived 
nebular conditions and abundances, of relevance to the 
present comparisons. 



slit nebular line intensities are taken directly from Esteban 
et al. (1991). Note that identical stellar ionizing flux dis- 
tributions to those obtained for the higher reddening case 
are used in this analysis, except that the absolute K-band 
magnitude is adjusted by 0.15 mag to Mk = —6.1 mag, 
for consistency with the adopted distance of 5 kpc. Re- 
sults from this alternative approach are very similar to 
those discussed here for Eb-v=1-3 mag. 

The electron density at the inner boundary was se- 
lected so that the predicted Ha flux agreed with observa- 
tions. For the Eb-v=0-9 mag case, typical quantities are 
iV e =l,350 cm~ 3 at the inner boundary (8" or 0.2 pc), and 
320 cm -3 at the outer boundary (50" or 1.25pc). Similar 
values were measured by Esteban et al. (1991), namely 
-1,000 cm" 3 at a distance of 10-25" and -200 cm" 3 for 
>30" (their region W). 



From Table |J, the discrepancy between the observed 
plasma conditions and photo-ionization modelling from 
the unblanketed CMFGEN analysis is more acute than Es- 
teban et al. (1993) obtained, for both reddening solutions. 
(A higher stellar temperature is imposed here by mod- 
elling the stellar spectrum, including the allowance for 
light metals, whereas Esteban et al. did not adopt any 
stellar constraint). 

In this case, the ionization balance in the nebula is 
much higher than observed, with S 2+ /S + overestimated 
by a factor of 40. In addition, the mean T e is far too high, 
with the predicted range from 11,200K at the inner nebula 
to 9,500K at the outer boundary, contrasting sharply with 
the observed value of 6,200K. Recalling Fig. a relatively 
strong flux at energies greater than the S 2+ edge at 356A 
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Table 4. Comparison between observed and predicted 
nebular properties of Ml-67 from CLOUDY, based on vari- 
ous ionizing flux distributions. Nebular line measurements 
are taken from slit position A of Esteban et al. (1991), de- 
reddened according to i?B-v=l-3 (Sect. [O] ) and given 
relative to H/3=100, including formal errors from Este- 
ban et al., neglecting systematic (calibration) uncertain- 
ties. I(H/3) represents the de-reddened, integrated H/3 flux 
received at the Earth (ergcm _2 s _1 ), obtained via the Ha 
flux measured by Grosdidier et al. (1998). 



Quantities 


Observed 


CMFGEN 


CMFGEN 


ISA-wind 


Blank? 




no 


yes 


yes 


log I(H/3) 


-9.7 


-9.7 


-9.7 


-9.7 


e 




0.05 


0.05 


0.05 


N e (r =0.2pc) 




2050 


1650 


2000 


iVe(r=1.25pc) 




480 


380 


460 


12+log N/H 


8.21 


8.21 


8.21 


8.21 


12+log O/H 


8.09 


8.09 


8.09 


8.09 


12+log S/H 


6.64 


6.64 


6.64 


6.64 


He/H 


0.2 


0.2 


0.2 


0.2 


3727 [0 ii] 


21±2 


37 


52 


52 


5007 [O in] 


<1: 


316 


0.0 


148 


5755 [Nil] 


0.7±0.1 


1.1 


1.9 


1.3 


5876 He I 


1.6±0.3 


30 


0.0 


31 


6548 [Nil] 


71±1 


24 


113 


45 


6563 Ha 


194±1* 


284 


292 


287 


6584 [N n] 


204±2 


71 


333 


134 


6717 [Sn] 


8.4±0.2 


1.8 


12.5 


2.7 


6731 [S ii] 


10.3±0.2 


1.9 


12.8 


2.8 


7135 [Arm] 


0.2: 


68 


0.1 


50 


9069 [S in] 


5.5±0.3 


30 


9.5 


24 


9546 P8 


1.9±0.2 


3.8 


3.9 


3.9 


N 2 +/N+ 




10 


0.00 


3 


o 2+ /o+ 


<0.02 


5.5 


0.00 


1.5 


s 2 +/s+ 


0.9 


35 


1.2 


17 


Ar 2 + /Ar+ 




57 


0.00 


25 


T.(K) 


6,200 


9,800 


6,900 


8,500 



*: See text for an explanation why the Ha flux does not agree 
with the Case B value using our stellar reddening. 



indicates that the predicted nebular ionization balance im- 
plies S 2+ » S 3+ > S + . Similarly for oxygen, the strong- 
flux predicted below the + edge at 353A implies an error 
in the 2+ /0+ ratio of order >300! Consequently, [Oni] 
A5007 is predicted to have an intensity that is greater than 
H/3, yet [O ill] A5007 is barely detected in Ml-67. Simi- 
lar problems relate to [Nil] A6584, Hei A5876 and [Arm] 
A7135. 

In contrast, the line blanketed CMFGEN model, with a 
much softer extreme UV energy distribution, produces a 
much lower predicted electron temperature, in the range 
6,600-7,900K, now in reasonable agreement with observa- 
tions. From Table ^, almost every nebular line is predicted 
to be within a factor of two of the de-reddened value, with 



the predicted ionization balance in good agreement with 
observations. 

An explanation for the success of this model can be 
understood by recalling Fig. ^| in which negligible flux is 
predicted at energies greater than the He edge at 504A, 
and indeed only a small flux beyond the S + A532 edge. 
Consequently, [S n] A6731 and especially [S in] A9069 are 
fairly well matched for either case, although the predicted 
ionization balance is actually marginally too low: S + ~ 
S 2+ ^ S 3+ . Another important success of this model is 
that negligible flux is predicted above the + 353A edge, 
so that the observed absence of [O in] A5007 is accounted 
for. 

However, the ionization of the model predicts negligi- 
ble He I A5876 and [Arm] A7135 emission, yet these are 
both observed, albeit weakly. Note that the precise choice 
of the He/H abundance ratio in the nebula does not play a 
major role in the predicted Hei A5876 strength, since it is 
primarily determined by the number of He I ionizing pho- 
tons. Less crucial deficiencies of this model are that [Nil] 
A6584 and [O n] A3727 are predicted to be too strong, al- 
though [Nil] /Ha (independent of reddening uncertainties) 
is in good agreement. 

5.2.2. IS A- wind 

Turning to ISA- wind, its energy distribution lies intermedi- 
ate between the unblanketed and blanketed CMFGEN mod- 
els, such that the nebular ionization structure is predicted 
to be somewhat higher than the observed properties, with 
T e overestimated by 40%. Once again, a significant ioniz- 
ing flux is predicted to lie at energies above the + edge, 
so that the predicted 2+ /0 + ratio is a factor of >75 
times too high, with very strong [Oni] A5007 predicted. 
Similarly, the flux above the S 2+ and He edges are too 
strong, so that that the S 2+ /S + ratio is in error by a factor 
of 20, and He I A5876 is strongly overestimated. 

5.3. Influence of uncertain nebular properties 

How dependent are the predicted nebular properties on 
our assumed filling factor, e, or geometry? As a test, we 
have carried out calculations for the blanketed CMFGEN 
case with £'b-v = 1-3 mag. 

For e=0.025 (a factor of two smaller), the required elec- 
tron density is \/2 higher, while predicted nebular proper- 
ties and line strengths properties are barely affected, with 
the electron temperature unchanged. The greatest influ- 
ence of a decrease in filling factor is a decrease in the 
sulphur ionization balance, from S 2+ /S + =1.2 to 0.9, with 
changes in nebular line strengths of up to 20% for [S H-Ill] . 
Comparable increases in S 2+ /S + and decreases in electron 
density are obtained for higher filling factors. 

Similarly, we have compared results with a constant, 
rather than exponentially decreasing, density through- 
out the nebula. For e=0.05, the required electron den- 
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sity is 7V e =530 cm -3 . Again, the influence is generally 
small, with the electron temperature unchanged and line 
strengths modified by up to 10%. Once again, the effect 
on sulphur is greatest, such that the ionization balance is 
decreased from S 2+ /S + =1.2 to 1.0 and changes in [Sll-Iil] 
strengths of up to 20%. 

From these calculations, we consider that the nebu- 
lar predictions are robust, for the blanketed CMFGEN case 
at least, except that the sulphur line strengths and ion- 
ization balance suffer from uncertainties of order ±50%. 
Equivalent tests for the unblanketed CMFGEN case indicate 
similar uncertainties for the ionization balance of oxygen 
and nitrogen, although the effect on emission lines from 
[O II — III] is minor. 

Finally, one further note of caution is necessary. 
Lamers et al. (1999) emphasise that nebular line fluxes 
and (absolute) abundances may be sensitive to clumping, 
of particular concern for the Ml-67 nebula. Although be- 
yond the scope of the present work, one would ideally like 
to carry out a re-evaluation of the nebular abundances and 
conditions of Ml-67, accounting for its highly clumped na- 
ture. 

5.4- Why are predicted ionizing fluxes different? 

A definitive explanation for the difference in ionizing spec- 
tra produced by the blanketed CMFGEN and ISA-wind 
models shortward of A700A is not straightforward, given 
the differences in their treatment of the radiative transfer, 
blanketing, atomic models and atomic data. Nevertheless, 
several differences in the alternative techniques deserve 
consideration. 

The most promising reason for the differences is that 
CMFGEN allows for photons that are absorbed in lines and 
re-emitted at other (longer) wavelengths. In contrast, no 
such channel is available in the is A- wind MC calculations, 
thereby overestimating the ionizing flux. Other possible 
explanations, such as photon loss operating in CMFGEN at 
the He I A584 resonance line, or an incorrect iron ionization 
balance being used for the MC calculation in ISA-wind, 
appear to be excluded (Sect. |3.2| ). 

In addition, the number of lines treated in both cases is 
incomplete, so that blocking will be underestimated, and 
the predicted ionizing fluxes overestimated. Recall that 
the number of lines considered by CMFGEN is relatively 
small, principally treating Fe in-v, while ISA- wind consid- 
ers many more lines and ionization stages, although this 
is clearly not complete (solely lines between measured en- 
ergy levels are included). Depending on the strength of 
this effect, softer ionizing energy distributions would re- 
sult, leading to improved agreement with nebular proper- 
ties for ISA-wind. 

In spite of the modelling and nebular difficulties, we 
conclude that the CMFGEN line blanketed ionizing en- 
ergy distribution for WR124 provides excellent consis- 
tency with the nebular properties of Ml-67, and solves 



the problem that was identified by Esteban et al. (1993) 
for this system based on unblanketed models. Use of ISA- 
wind leads to consistent stellar properties for WR124, al- 
though the consistency with observed plasma conditions 
for Ml-67 are somewhat poorer than for CMFGEN. 

6. Conclusions 

We have used WR124 and its associated nebula Ml-67 to 
assess the reliability of Lyman ionizing flux distributions 
from current stellar models. For WR124, we find that stel- 
lar properties obtained with two line blanketed, non-LTE 
extended model atmospheres, ISA-wind (de Koter et al. 
1997) and CMFGEN (Hillier & Miller 1998) are in excellent 
agreement. However, comparisons with the observed prop- 
erties of Ml-67 favour the ionizing energy distribution ob- 
tained by CMFGEN rather than ISA-wind, reconciling pre- 
vious inconsistencies identified by Esteban et al. (1993). 
Certainly, more general conclusions await comparisons for 
a variety of stars. 

In the future we propose to extend our study to WR 
stars of different spectral type. One well known exam- 
ple is HD 192163 (WR136) and its ring nebula NGC 6888, 
for which Esteban et al. (1993) also identified a major 
inconsistency. However, in contrast with Ml-67, there is 
evidence that nebular emission in NGC 6888 may be pro- 
duced by shocks rather than photo-ionization. Further, 
the properties of other Galactic WR nebulae are often 
poorly known, such that different stellar ionizing distribu- 
tions can reproduce the nebular properties by selecting, 
for example, alternative electron densities or shell thick- 
ness. WR nebulae in the LMC may provide more suitable 
candidates (Dopita et al. 1994). 

In general, study of young stellar populations in ex- 
ternal galaxies via associated H II regions relies heavily on 
suitable stellar ionizing flux distributions (e.g. Lcithcrer et 
al. 1999). Currently, the only WR models that are gener- 
ally available at this time are unblanketed, pure helium 
energy distributions (Schmutz et al. 1992). The calcu- 
lation of a large multi-parameter grid of line-blanketed 
models remains a formidable computational challenge, un- 
less codes such as ISA-wind can be relied upon. Crowther 
(1999) finds that, in general, the Hi and Hei ionizing 
fluxes of early-type WR stars are not greatly affected by 
line blanketing. However, the additional blanketing from 
C and O in early WC stars produces a softer ionizing spec- 
trum than for early WN stars of identical temperatures, 
with negligible flux emitted A <300 A. Indeed, unblan- 
keted models are probably to be preferred for the study 
of young starburst regions in very metal deficient regions. 
However, the (uncertain) role that metallicity may play 
on mass-loss rates of Wolf-Rayet stars does need to be 
quantitatively investigated (e.g. Crowther et al. 1999). 
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